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The Callan-Treiman low-energy theorem offers an opportunity to test electroweak couplings of 
light quarks to the gauge boson W. To that aim, we introduce a model-independent and accurate 
dispersive parametrization of the two form factors. We then discuss three applications to 
the analysis of K^t, and K^^t, measurements: the prediction of the ratios Y {^K^^^IYiKg-},), the ex- 
traction of I/+ (0)y„.s| and finally the possible measurement of m„ — m^i induced isospin breaking 
asymmetry. 
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1. Introduction 

The hadronic matrix element involved in K(/i, decays is described in terms of two form factors 
/+ and The vector form factor f^^{t) represents the P-wave projection of the crossed channel 
matrix element {0\s'Y^u\Kn) , whereas the S-wave projection is described by the scalar form factor 

/i'"(0 = /f + , (1-1) 

mi -mi 



where t = [pji — Pk)^. In the following, we will consider the normalized form factors /o(f) 
/f ''"(O), f+it) = /f ^"(0//f ''"(0), /o(0) = /+(0) = 1, and we will try to descril 

their shape as precisely as possible in the region of interest, in order to test the Standard Model. 



2. Callan-Treiman theorem: a test of the Standard Model 

The low energy theorem of Callan and Treiman predicts the value of /o(0 at the Callan- 
Treiman (CT) point, namely t = Akti = mj^ — m|, in the SU(2) x SU(2) chiral limit (mu,md 0). 
When mu,md ^ 0, one writes 

C^MA..) = %j^^^^, (2.1) 

where the CT correction Act ~ ^ {mu.d/4nFji). This correction has been estimated within Chiral 
Perturbation Theory (ChPT) at next to leading order (NLO) in the isospin limit [Q] with the result 

A^^^ = -3.5 10"^ . (2.2) 

Note that the CT point Akti is situated between the physical end point of K^i, decays to = {rriK — m^j)^ 
and the Kn scattering crossed channel threshold Ikk = {mK + mjiY . Fk^k are respectively the kaon 
and pion decay constants parametrizing QCD effects. Act is not enhanced by any chiral logarithm, 
so its value is rather small. A complete study of ff{p^) corrections [Q] is not available yet. Howe- 
ver, we expect that these effects should not change the value of A^lfP in Eq. ( ^!^ by one order 
of magnitude. Concerning the study of isospin breaking corrections, contributions from the tt^-tj 
mixing in the final state, usually responsible for large isospin violations, are absent in the neutral 
K^-i mode, K*^ ^ . In the charged mode instead, these effects can easily reach a few percents. 
Furthermore, the electromagnetic (EM) effects have still to be fully investigated. 

Summarising, Eq. ( p] ) allows to predict with good precision the value of /()(?) at the CT point 
and test the Standard Model (SM) couplings in a sector where they have not been tested so far. 



More specifically, assuming the SM couplings, we can test Eq. (gjj) by measuring independently 



the l.h.s and r.h.s of Eq. {lA): first step, in the SM the QCD parameters on the r.h.s of Eq. ( |2. 1[ ) can 



be extracted from experimental information, such as the ratio r^+^ 101' \ f+^ (0)V„.sI [| 

from K^e^Ve and Vud [Pi, namely 



CsM ■ 



Fk^Vus 



ud 



FK^n^(C^\u ' \+^CT= Bexp + Act , (2.3) 



\fr{o)Vus\ 
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with Bexp = 1 .2438 ± 0.0040 . It can be noticed that this prediction is independent of the knowledge 



of Vus- Besides, since experimentally one measures \FK+yus/pK+ Vud \ , this will introduce in Eq. ( p.l[ ) 
some small isospin breaking effects due to the difference between Ff^+ and Finally, the l.h.s 



of Eq. ( |2.1[ ), f(){m\ — m\), can be determined with some care from the K^j, decay distribution. 



To this end, one has to extrapolate /o from the physical region (O < f < ?o = {mK — niji)'^) up 
to Akk- Fits to the K^^ measured distributions are usually performed using a Taylor expansion of 
the two form factors 



/o,+ (0 = l+A(: 



:>,+ 



mi, 



+ ... 



(2.4) 



So far, the results available for the scalar form factor only give the slope X^^'^, since the different 
experiments cannot discriminate Ag^^ from the curvature which are highly correlated. The results 
disagree among all the experiments as illustrated in Fig.[l|a. To extrapolate up to Akti, we cannot 
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Figure 1: a) Shape of /o with a linear fit (Aq^'') [0 compared to the dispersive analysis, b) Icr domains 
allowed by the branching ratio measurements and by the NA48/DR analysis. 

neglect the positive curvature. Anyway, even if there was no curvature, the SM prediction for the 
slope implied by the CT theorem would be Ao = 0.0208 ± 0.0003 and that is above all the exis- 
ting experimental results. Note that the ChPT prediction for Ao is much less precise than the CT 
theorem since it derives from an expansion in powers of m^. However, the curvature being positive, 
it is obvious that the measured slope X^^'^ can be at most an upper bound for the "true" slope 
mathematically defined as: Ao = m|/'(0). Under these conditions, X^^'' may depend on the fitted 
distributions and thus on different weights given to the different areas of the Dalitz plot by various 
experiments. That could explain the apparent discrepancy between the experimental results. 

3. Dispersive representations of the Kn scalar and vector form factors 



3.1 Scalar form factor 

Our aim is to construct a very precise representation of the scalar form factor between and 
Afcjt in order to test the SM prediction at the CT point. The information we have is: the value of 
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(3.1) 



/o at zero, the Ktt scattering phase in the isospin limit and the fact that /o has to vanish as tff{l/s) 
for large negative s [|n|]. In order to determine C, one can write a dispersion relation for ln(/o(f)) 
twice subtracted at and 

"^(lnC-G(0)l,withG(0 = ^^^^^^^ P , 

-^Kn -I JtKn S {s-AKK){s-t-ie) 

where tKK is the threshold of Kn scattering and 0(s) is the phase of fo{s). We decompose the 
integration range into two parts: the elastic part (Iktc < s < A) and the (inelastic) asymptotic part 
(A < 5 < oo) corresponding to the splitting G{t) = GK7t{A,t) +Gas{A,t) ± 5G{t) , where dG repre- 
sents the total uncertainty. For the elastic part, GKn{^,t), Watson theorem states that in the elastic 
region and in the isospin limit (j){s) equals the S-wave, 7=1/2 Ktt scattering phase 5k7i{s). This 
phase has been extracted (in the isospin limit) from experimental data [|T^] solving the Roy-Steiner 
equations [13]. For the asymptotic part, Gas{A,t), we take 0(5) = ^as{s) = 71 ±71 since 0(+oo) = 71 
due to the behaviour of /o at large negative s. Note that we consider a very conservative esti- 
mate of the asymptotic uncertainty. Thanks to the two subtractions, G{t) in the region of interest 
is relatively insensitive to the unknown phase at high energy and also to the precise knowledge 
of A, characterizing the end of the elastic region. In Fig.|2| is represented the function G{t) for 
A = 2.77 GeV^. Let us stress that the uncertainties on G{t) represent at most 10% of its central 
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Figure 2: G(f) with the uncertainties dGcs and 
5Gkji added in quadrature. 
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Figure 3: H(t) with the uncertainties 5Has and 
5Hkk added in quadrature. 



value (G(0 < G(0) = 0.0398 ± 0.0040) which by itself does not exceed 20% of the expected value 
of InC. Hence, one can be confident that the theoretical uncertainty is under control and will not 
be the dominant error in the experimental extraction of InC. We have sticked here to the isospin 
limit. The potentially dangerous contribution of the cusp at (m^+ +mj^o)'^ due to isospin breaking 
should however not affect our results. Besides, a zero of /o(f) could only influence our dispersive 
construction if it was located at t <tKn, close to the physical region. This seems unhkely due to 
/o(0) = 1 and the smallness of the slope. For more details on these two last points see Ref. [^. 
In order to determine InC, in principle there exists a sum rule 



G( 



=) = InC 



(3.2) 



dictated by the asymptotic behaviour of /q. However, it is not precise enough to allow to determine 
InC with the needed accuracy without adding any information on the high energy behaviour of 
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the phase of the form factor. Indeed, the integrand defined in Eq. (3.1) drops as 1/^^^^ whereas 
the one in G{—°°), which has one less subtraction, drops as \/s^. For instance, if we vary A 
between 2.25 GeV^ and 2.77 GeV^, G(-oo) varies between 0.1335 and 0.3425, whereas G{t) 
is practically not affected. We will thus instead determine InC from experiment by fitting the K^j, 
decay distribution with the dispersive representation formula of /o(f), Eq. (^?T|). Once InC is known, 
we can estimate the first two coefficients of the Taylor expansion, Eq. (^^, and their correlation 

2 4 

Ao = ^(lnC-G(0)) and = Aq^ _ 2^G'(0) = Aq^ + (4.16 ±0.50) x 10"^ (3.3) 

3.2 Vector form factor 

One can use the same dispersive construction as for the scalar form factor. Writing a twice 
subtracted representation of the form factor but this time at (/+(0) = 1 and /+(0) = A+/m|) 
leads to 

A(,)=exp[i,(A, + H(,))] . where Hi,) £^4,^^^. (34, 



Here, in the elastic region, the phase of the vector form factor 9(5) equals the 7 = 1/2, P-wave 
scattering phase. The procedure using Roy-Steiner equations suffers in this case from a lack of 
relevant experimental inputs. Hence, we have constructed a partial wave amplitude using a Breit- 
Wigner parametrization around the A"* (892) pole (using the PDG value as input) which is analytic, 
unitarized and has the correct threshold behavior following Ref. [|T4{]. The resulting function H 
is decomposed into two parts as previously with A ~ (1.4 GeV)^, see Figj3[ Note that in this 
channel the A'* (892) pole dominates and the vector form factor is very well described by the pole 
parametrization, which is in full agreement with our dispersive construction. For other works on 



the subject, see Ref. [15] and references therein 



3.3 First dispersive analysis 

NA48 has realized the first dispersive analysis (referred as NA48/DR in the following) of their 
/:^3 Dahtz distribution leading to the result InC'^'P = 0.1438(140) and A'^'' = 0.0233(9). 
The sum rule //(— 00) = — A_|_ is well fulfilled. However, in the case of InC, InC^^ is close to the 
lower bound of G(— 00). Even if the sum rules are not very stringent since they are sensitive to 



the behaviour of the phase at high energy, the sum rule, Eq. (3.2), indicates that the phase of /o(0 
should drop after A before reaching its asymptotic value 7i. Note that a similar phenomenon exists 
in the case of nn scattering [17]. The result for /o(0 is shown in Fig.|I|a for the SM prediction. 



Eqs. (gj) and ( |2.2| ), and for the NA48/DR result. The NA48/DR analysis indicates a 5a deviation 
with the SM prediction that could be interpreted as a signal of physics beyond the SM as for 
instance a first evidence of a direct coupling of W to right-handed quarks [|], [I^. However, the 
apparent disagreement between the different experimental measurements has to be clarified before 
drawing a firm conclusion. 



According to Eq. (3.3), one can infer in the case of the SM value, Eq. (2.3), Aq"^^ = 0.01523 ± 



0.00046 +0.069Act^ A^™ = (6.48 ±0.52 + 47. lAcr) 10""^ and in the case of the NA48/DR 
result : Aq^^^s ^ o.0089 ± 0.0012, A^ ^'^^s ^ (4 95 ^ o.55) 10"^ . 



This is in complete agreement with the prediction of [19p 



5 



Dispersive representation and shape of Kfi, form factors 



Emilie Passemar 











l/+(0)K.lf 


i/+(o)K.v!r 


a'' 


SM 


0.6677(28) 


0.6640(28) 


SM 


0.21642(60) 


0.22176(100) 


2.47(54)% 


NA48/DR 


0.6589(33) 


0.6552(33) 


NA48/DR 


0.21834(69) 


0.22378(107) 


2.49(54)% 



Table 1: Left: Prediction of R^j^ ■ Right: Extraction of |/+ (0)V„.v | and Asij^2) for the muonic mode using 
calculated with the dispersive approach together with the experimental inputs from [ p^ ] and the preliminary 
updated EM corrections from [pO[|. 



4. Applications 

The knowledge of /o(f) and f+{t) allows to determine the phase space integrals /^^^q entering 
the well-known master formula of the Kij partial widths, intensively discussed during the confe- 
rence (see Talks in the Vus session), 

r^+/o = ^+/o Sew (1 + 2A™/o^) \f+^'\o)Vus\^ 4^/o , A+/o = 4^/0 Cj m|^/o/(1927r3). (4.1) 

Here, we do not consider the isospin breaking (IB) correction, A5j/(2) = /+ ''^''(0)//:^*''^ (0) — 1, 
mostly induced by m„ — effects, as a known theoretical input since we would like to determine it 
experimentally. Assuming jj^/e universality, separately tested in the ratio r(^e2)/r(^^2) [H|], the 
knowledge of can be converted into a prediction for R^j^ = ^K+zOjiif^K+zoeS- Using the last 
update of the EM corrections [^], we obtain the results in Tab. 1 (left). For the evaluation of /^+/o, 
we have used in the case of the SM prediction, InC from Eq. (2.3) with Act from Eq. (2.2) and 
A+ = 0.0245, corresponding to the K*{S92) pole, and in the case of the NA48/DR result, InC^^ 
and A"^^ including the correlation from this analysis. These calculations can be compared with 
the experimental results summarized in the legend of Fig. lb. R^j'^ depends almost exclusively on 
InC, since the dependence on A+ partially cancels by taking the ratio. This offers the possibility of 
extracting InC quasi independently of A+ contrary to the Dalitz plot analysis. In Fig.|l]b are plotted, 
in the plane (InC, A+), the la domains allowed by the individual /T" and branching ratio (BR) 
measurements. Note that one loop ChPT calculations [^] show that the form factor shapes (InC, 
A+) are, to a good approximation, not affected by IB. For comparison, we show on the same plot 
the independent result of the NA48/DR analysis, which is not incompatible with the BR measure- 
ments, though on the lower side of the domain. 

Using the master formula, Eq. (p~T|), we can extract |/+(0)V'„i| from the experimental mea- 
surements (see the Flavianet fit [|I0|]). Waiting for other experiments to perform a similar dis- 
persive analysis, the extraction from the muonic mode, Tab.l(right), only involves the NA48/DR 
result and hence has still a limited accuracy. As it has been emphasized during the conference, see 



Ref. [|22|], the measurements of [/+(0)V'„.v| from the neutral and charged modes as well as the EM 



corrections |gOp have now reached a sufficient degree of accuracy to give an opportunity to test the 
theoretical prediction A5(/(2) = 0.0231(22) [ p3| ] and the input value of R = {nis — rh)/ (mj — m„). 
As it has been shown in Ref. [|^] and from Tab. 1 (right), using an average of the BR measure- 
ments, A5[/(2) riow seems closer to the theoretical prediction. 



In conclusion, we have introduced a physically motivated model-independent dispersive re- 
presentation of the Kf^i, form factors which only involves one free parameter. This will improve 
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the accuracy of the decay analysis and can help to partially remove the apparent discrepancy 



between the different experimental results (see also the discussion in Ref. [g4j]). It will then allow 
to test the SM and to extract |/+(0)V„i| with a better precision. 
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